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SIGNAL  SIGNATURES  OF  TOPOGRAPHIC 
FEATURES  USING  ANALOG  TECHNOLOGY 


INTRODUCTION 


The  objective  of  this  research  was  to  develop  signal  signatures  and  signatures 
ot  spectral  decompositions  of  selected  man-made  topographic  features,  such  as  linear 
roads,  road  intersections,  and  rectangular  buildings.  An  analog  signal  processor  was  pro¬ 
posed  as  the  signature  detection  device  because  of  its  high  speed  and  large  bandwidth 
capability. 


The  gray  shade  distribution  of  selected  topographic  features  from  aerial  photographs 
was  converted  into  analog  electronic  signals  by  an  area  sensor  array.  The  signals  were 
passed  through  the  threshold  gate  to  eliminate  unwanted  background  noise  and  were  con¬ 
structed  as  two-dimensional  binary  pictures  that  represent  spatial  signal  signatures  of  the 
selected  topographic  features.  The  threshold  signals  were  transformed  into  the  Walsh 
domain,  and  various  spectral  components  were  obtained.  An  electronic  system  capable 
of  generating  two-dimensional,  32  x  32  low-order  Walsh  funtions  for  decomposing  topo¬ 
graphic  feature  images  was  designed  and  built. 


The  signal  signatures  were  compared  in  both  domains  to  identify  uniqueness,  sim¬ 
plicity,  and  detectability.  A  scheme  for  detecting  the  decomposed  spectral  components 
of  the  selected  topographic  features  using  an  analog  processor  is  presented.  In  the  last 
sections,  conclusions  based  on  the  research  are  listed. 


PRESENT  SYSTEM  CONFIGURATION 


A  system  that  will  produce  signal  signatures  from  the  gray  shade  distribution  of  a 
photo  transparency  in  a  short  time  and  with  recognizable  resolution  was  conceived  and 
developed.  The  block  diagram  of  such  a  system  is  shown  in  figure  1. 
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Figure  1.  System  Block  Diagram. 


A  9-  by  9-inch  aerial  transparency  is  illuminated  by  a  white  light  source,  and  a  section 
of  the  image  is  projected  on  a  Reticon  32-  by  32-element  area  sensor  array  through  a 
lens.  The  array  converts  the  optical  energy  of  the  image  into  a  video  signal.  The  video 
signal  going  through  the  interface  electronic  is  passed  through  the  threshold  gate  of  a 
Hewlett-Packard  (HP)  2100  mini-computer  and  becomes  two  values,  “on”  and  “off.”  The 
value  of  the  threshold  is  variable,  which  provides  a  very  convenient  means  for  isolating 
signals  representing  the  selected  feature  image  from  the  unwanted  background  noise.  The 
“on”  value  of  the  video  is  set  to  “100"  and  “off  to  “0,"  and  the  set  values  are  printed 
by  an  HP  line  printer  as  a  two-dimensional  binary  array  of  32  x  32  pixels  that  represent 
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the  spatial  signal  signatures  of  the  selected  topographic  features.  After  passing  through 
the  threshold  gate,  the  signal  is  also  stored  in  32  x  32  memory  locations.  This  stored 
signal  is  multiplied  with  a  two-dimensional  Walsh  function  that  has  two  values,  either 
+  1  or  -1  in  magnitude,  and  it  has  the  same  32  x  32  clock  periods  as  the  input  signal. 
The  multiplied  outputs  for  each  frame  are  added  and  the  total  becomes  the  Walsh  trans¬ 
form  coefficient  for  that  particular  Walsh  function  used.  This  process  implements  the 
two-dimensional  discrete  Walsh  transform,1 


32  32 

a (i, j)  =  — - k  =  l  1  =  1  f(k,  1)  Wal(i.k)  Wal(j,  1) 

1024 


where  f  (k,  1)  is  the  image  binary  array,  Wal  (i,  k)  Wal  (j,  1)  is  the  two-dimensional 
Walsh  function  of  order  i  and  j,  and  k,  1,  i,  and  j  take  values  from  1  through  32.  By 
sequentially  changing  the  order  of  the  Walsh  function,  the  complete  set  (32  x  32)  of  the 
Walsh  transform  coefficients  was  obtained.  The  walsh  functions  used  were  generated  by 
the  minicomputer.  Figure  2  shows  the  first  8 -by  8-order,  two-dimensional  Walsh  fun¬ 
ctions  in  the  range  of -Vi  to  +Vz.  The  black  areas  indicate  +1  in  magnitude, and  white -1. 


Since  Walsh  transform  coefficients  were  produced  by  using  Walsh  functions  having 
spatially  alternating  magnitudes  (from  +1  to  -1)  at  different  sequences,  the  Walsh  trans¬ 
form  coefficients  are  also  decomposed  spectral  components  of  the  signal  signature  of  the 
input  signal.  These  spectral  components  were  divided  by  a  convenient  constant,  in  this 
case  1 .024,  which  is  the  number  of  pixels  in  a  frame  for  normalization  purposes.  The 
entire  normalized  set  of  the  spectral  components  was  printed  out  by  the  line' printer. 


'h.  p.  Harmuth,  Sequency  Theory,  Foundation  and  A  p  plication ,  New  Y  o  rk  ,  A  cadem  ic 
Press,  Inc.,  1977,  pp.  55-5  6. 
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ELECTRONIC  WALSH  FUNCTION  GENERATOR 


The  signal  signatures  of  the  decomposed  spectral  components  described  in  the  pre¬ 
vious  section  were  obtained  by  using  the  minicomputer  with  the  video  signals  from  the 
array  as  input.  It  takes  approximately  2  seconds  to  produce  a  complete  set  of  32  x  32 
Walsh  transform  coefficients.  To  reduce  this  long  processing  time,  a  hardware  electronic- 
system  capable  of  generating  two-dimensional,  32  x  32,  lower  order,  Walsh  functions  was 
designed  and  built. 

Figure  3  shows  the  block  diagram  of  this  Walsh  function  generator.  The  entire  gener¬ 
ator  consists  of  seven  subfunctions  as  shown.  Since  the  X  function  generator  and  the  X 
decoder  are  functionally  identical  to  the  Y  function  generator  and  the  Y  decoder,  only 
one  pair  of  them  will  be  discussed.  The  X  function  generator  consists  of  a  5-digit  binary 
counter,  26  exclusive  OR  gates,  and  6  inverters  as  shown  in  figure  4.  The  binary  counter 
(a  SN7493  4-digit  binary  counter,  plus  a  SN74H106  flip-flop)  uses  a  modified  clock  as 
input  to  generate  a  5-digit  binary  counting  sequence  A,  B,  C,  D,  and  E.  By  properly 
gating  one  of  these  digits  with  another  or  with  an  already  produced  lower  order  Walsh 
function,  the  first  32  lower  order,  X-direction  Walsh  functions,  Wal  (0,  x)  through  Wal 
(31,  x)  are  generated.  The  end-of-lines  (EOL)  pulse  from  the  array  resets  the  whole  se¬ 
quence  after  32  clock  pulses,  and  the  next  following  modified  clock  restarts  the  whole 
sequence  again.  This  process  repeats  until  a  set  of  decomposed  signal  signatures  is  ob¬ 
tained. 


The  Y  function  generator  is  functionally  identical  to  the  X  function  generator  except 
that  the  input  of  the  binary  counter  is  the  EOL  pulses  and  the  reset  pulses  are  the  end- 
of-frame  (EOF)  pulses.  This  generator  produces  Wal  (0,y)  through  Wal  (3 1  ,y)  as  shown  in 
figure  5. 


The  X  decoder  consists  of  a  5-digit  binary  counter,  a  dual  2-to-4-line  demultiplexer, 
four  dual  4-to-l-line  multiplexers,  a  reset  monostable  multivibrator,  and  eight  OR  gates 
(see  figure  6).  The  binary  counter  (a  4-digit  binary  counter  SN7493,  plus  a  flip-flop 
SN74H106)  uses  the  frame  reset  pulses  as  input  to  generate  a  5-digit  counting  sequence 
A,  B,  C,  D,  and  E.  Digits  C,  D,  and  E  are  used  as  input  to  the  dual  2-to-4-line  demulti¬ 
plexer  for  generating  Gj  to  Gg  pulses.  The  main  decoding  task  is  done  by  four  dual 
4-to-l-line  multiplexers  (SN74153).  The  X  Walsh  functions  generated  by  the  X  function 
generator  are  permanently  connected  to  32  inputs  of  these  multiplexers,  as  shown  in 
figure  6.  By  proper  selection  of  A,  B,  and  G,  through  Gg,  Wal  (0,x),  Wal  (1  ,x),  Wal  (31  ,x) 
will  appear  sequentially  in  one  of  the  outputs  of  the  multiplexers  for  each  frame  period 
of  the  array.  The  outputs  of  the  multiplexers  1  to  8  are  gated  through  eight  OR  gates  to 
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Figure  3.  Block  Diagram  of  a  Walsh  Function  Generator. 
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Figure  6.  X  Decoder. 
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become  a  logical  one-line  output  of  Wal  (i,x)  as  shown.  At  the  end  of  decoding  the  32  X- 
direction  Walsh  functions,  the  binary  counter  generates  a  Yjn  pulse  to  be  used  as  input 
for  the  Y  decoder.  This  Y;„  pulse  also  resets  the  binary  counter  of  the  X  decoder.  The 
decoder  can  also  be  reset  to  its  starting  condition  anytime  a  reset  pulse  from  the  control 
circuit  occurs. 


The  Y  decoder,  shown  in  figure  7,  is  functionally  identical  to  the  X  decoder  except 
the  input  to  the  binary  counter  is  Yin  pulses  and  the  pulses  are  self  resetting  pulses. 
The  decoder  can  also  be  reset  to  its  starting  condition  anytime  a  reset  pulses  from  the 
start  circuit  occurs.  This  decoder  multiplexes  Wal  (0,y)  through  Wal  (0,31)  into  a  single 
line  of  Wal  (j,y). 


The  start  circuit  for  the  Y  decoder  consists  of  two  flip-flops  (SN74H106),  five  NAND 
gates,  an  NOR  gate,  an  inverter,  and  a  switch  (see  figure  8).  The  circuit  takes  in  deleted 
clock  pulses  as  input  and  produces  the  control  clock  as  shown.  The  control  clock  is 
produced  only  with  the  consecutive  appearance  of  both  the  frame  reset  from  the  array 
and  the  pressing  of  the  temporary  switch  S.  Once  started,  the  circuitry  will  continue  to 
produce  control  clock  pulses  until  all  32  Y-direction  Walsh  functions  are  decoded  and  the 
A'  pulse  generated.  This(A)  pulse  produces  the  reset  pulse  in  the  reset  circuit,  which  in 
turn  resets  the  start  circuit. 


The  clock  divider  consists  of  a  4-bit  binary  counter  SN7493,  a  flip-flop  SN74H106, 
three  NAND  gates,  a  3-input  NAND  gate,  and  two  switches  (see  figure  9.)  By  appropriate 
position  selection  ot  Sj  and  S2,  the  incoming  control  clock  pulses  will  be  divided  by  1 , 2, 
4,  and  8  as  desired.  The  divided  pulses  are  designated  as  modified  pulses.  The  reset  pulse 
resets  the  divider  to  its  initial  condition  for  new  incoming  pulses. 


The  reset  circuit  consists  of  two  portions  (figure  10).  In  the  first  portion,  the  reset 
circuitry  includes  two  NAND  gates,  an  NOR  gate,  an  inventor,  and  a  switch  S3  to  pro¬ 
duce  reset  and  reset  pulses  whenever  the  temporary  switch  S3  is  pressed.  In  the  second 
portion,  two  inverters  invert  the  Frame  Reset  and  EOL  (end-of-line)  pulses  from  the 
array  into  the  Frame  Reset  and  EOL  as  shown  in  figure  10. 


The  Walsh  function  generator  produces  X-direction  and  Y-direction  Walsh  functions 
on  a  pair  of  separate  lines.  In  the  following  section,  a  scheme  is  discussed  that  will  use 
outputs  of  the  Walsh  function  generators  as  two-dimensional  Walsh  functions  in  conjunc¬ 
tion  with  a  sensor  array  system  for  producing  Walsh  transforms.  The  Walsh  function 
generator  discussed  above  was  designated,  built,  and  tested  to  verify  feasibility  of  per¬ 
formance. 


Figure  9.  Clock  Divider. 
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FUTURE  SYSTEM  CONFIGURATION 


The  future  system  will  use  the  hardware  Walsh  function  generator  in  conjunction 
with  a  32-by  32-element  sensor  array  system  (see  figure  1 1.)  The  Wal  (i,x)  and  Wal(j,y) 
from  the  Walsh  Function  Generator  will  be  passed  through  an  exclusive  OR  gate  to  pro¬ 
duce  pulses  representing  +1  values  of  the  two-dimensional  Wal  (i,x)  Wal  (j,y ).  This  signal 
is  inverted  to  yield  -1  values  of  Wal  (i,x)  Wal  (j-Y )-  The  topographic  images  sensed  by  the 
array  (video  signal)  will  pass  through  the  threshold  gate,  MC14I4L,  as  shown.  The  signal, 
after  passing  through  the  threshold  gate,  is  gated  with  signals  representing  +1  and  -1  of 
Wal  (i,x)  Wal  (j,y)  in  separate  NAND  gates  with  the  inverted  modified  clock  as  controlling 
signals.  The  output  from  each  of  these  NAND  gates  is  counted  by  a  counter  for  an  array 
frame.  The  difference  between  these  two  numbers  (positive  or  negative)  represents  the 
Walsh  transform  coefficient  of  the  image  for  that  particular  Walsh  function.  The  out¬ 
puts  of  the  counters  are  converted  to  analog  signals  by  two  digital-to-analog  converters. 
The  converted  analog  signals  are  subtracted  in  a  differential  amplifier.  The  output  of  the 
differential  amplifier  represents  a  Walsh  transform  coefficient  of  the  topographic  feature 
sensed  by  the  sensor  array  for  a  Walsh  function.  By  sequentially  changing  i  and  j  of 
Wal  (i,x)  Wal  (j,y),  all  32  x  32  Walsh  transform  coefficients  will  be  obtained.  By  using 
the  hardware  Walsh  function  generator,  the  time  required  for  processing  Walsh  transforms 
is  much  less  than  the  time  required  by  using  software. 


A  SUGGESTED  DETECTION  SCHEME 


In  this  section,  a  scheme  for  detecting  the  decomposed  spectral  components  (or  Walsh 
transform  coefficients)  using  an  analog  processor  is  described.  In  the  scheme,  two  cas¬ 
caded,  32-stage  programable,  binary-analog  correlators2  are  used  (figure  12).  In  the  next 
section,  it  will  be  shown  that  the  significant  spectral  components,  in  most  cases,  appear 
in  the  first  8  by  8,  lower  order,  Walsh  transform  coefficients.  Thus,  two  cascaded.  32- 
state,  binary-analog  correlators  should  be  sufficient  to  detect  the  important  coefficient 
to  yield  recognizable  results.  The  Walsh  transform  coefficients  will  be  connected  to  the 
input  of  the  analog  delay  line,  and  the  binary  pattern  that  represents  known  features 


2V.  Strusilla,  A  Program  mahle  B  inary  -A  nalog  Correlator,  Rcticon  Corporation. 
Sunnyvale,  CA,  Technical  Note  No.  106 
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will  be  connected  to  the  input  of  the  static  shift  register.  Each  stage  of  the  analog  delay 
line  has  a  pair  of  taps.  These  taps  have  switches  in  series  with  them  that  are  controlled 
by  the  true  and  complement  outputs  of  the  static  digital  shift  register.  By  loading  a 
binary  work  into  it,  the  static  shift  register  will  select  the  taps,  which  are  connected  to 
two  output  lines,  thus  proving  the  ability  to  do  correlation.  By  sweeping  the  known 
binary  words  with  megaHertz  rate,  the  Walsh  coefficients  representing  a  particular  type 
of  feature  will  be  quickly  detected  in  the  output  of  the  operational  amplifier. 


COMPARISON  OF  SIGNAL  SIGNATURES  IN  BOTH  DOMAINS 


As  the  preliminary  test,  a  number  of  selected  manmade  features,  such  as  road  inter¬ 
sections  and  portions  of  straight-line  roads  from  aerial  photographs,  and  a  rectangular 
mask  pattern  were  used  as  inputs  for  the  present  system.  The  spatial  signal  signature  and 
the  associated  spectrally  decomposed  components  (Walsh  transform  coefficients)  for 
these  test  features  are  shown  in  figures  13  through  15.  It  is  seen  that  the  signal  signature 
of  the  spectrally  decomposed  components  is  much  simpler  than  the  spatial  signal  sig¬ 
nature  of  features  for  all  the  cases  shown.  In  most  of  the  cases,  the  significant  spectral 
components  are  distributed  among  the  few  lower  order  Walsh  transform  coefficients. 
Further,  each  transform  pattern  is  unique,  and  it  can  be  easily  distinguished  from  the 
rest.  For  example,  the  road  intersection  resulted  in  large  first-row  and  first-column  co¬ 
efficients,  and  the  horizontal  straight-line  road  yielded  large  coefficients  in  the  first 
column.  However,  the  rectangular  image  is  transformed  into  eight  significant  coefficients 
intersecting  rows  1  and  3,  and  columns  1,  3,  5,  and  7.  If  the  straight-line  road  is  turned 
90  degrees,  then  the  corresponding  significant  spectral  components  will  appear  in  the 
first  row  as  shown  in  figure  16.  The  signal  signatures  of  the  three  types  of  topographic 
features  considered  above,  after  turning  45  degrees  with  respect  to  x-y  coordinates, 
are  shown  in  figures  17  to  20.  Again,  a  very  simple,  yet  unique,  distribution  of  spectral 
decomposition  was  observed  for  all  cases.  Evidently,  a  simple  decision  procedure,  either 
in  hardware  or  software,  can  be  implemented  to  detect  and  recognize  the  above  limited 
set  of  features  automatically. 
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Figure  13.  Spatial  and  Decomposed  Signal  Signatures  of  a  Road  Intersection. 
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Figure  16.  Spatial  and  Decomposed  Signal  Signatures 
of  a  Vertical  Line  Road. 
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Figure  18.  Spatial  and  Decomposed  Signatures 
of  a  Diagonally  Oriented  Line  Road  ( 1 ). 
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Figure  19.  Spatial  and  Decomposed  Signal  Signatures 
of  a  Diagonally  Oriented  Line  Road  (2). 


Figure  20.  Spatial  and  Decomposed  Signal  Signatures 
of  a  Diagonally  Oriented  Rectangular  Object 


CONCLUSIONS 


1 .  The  signal  signature  of  the  spectrally  decomposed  topographic  features  is  much 
simpler  in  distribution  compared  to  the  spatial  signal  signature  of  the  same  topographic 
feature  for  all  selected  cases. 

2.  In  most  cases,  the  sigr.ift  roiif  spectral  components  are  distributed  among  few  lower 
order  Walsh  transform  coefficients  Further,  each  pattern  transform  is  unique  in  itself, 
and  it  can  be  easily  distinguished  from  the  rest. 


3.  A  simple  decision  procedure,  either  in  hardware  or  software,  can  be  implemented 
to  detect  and  recognize  the  selected  set  of  topographic  features  automatically. 


4.  The  present  sensor  array  and  minicomputer  system  requires  about  2  seconds  to 
obtain  32  x  32,  lower  order,  Walsh  transform  coefficients.  With  the  in-house  designed 
and  built  hardware  Walsh  function  generator,  this  process  time  can  be  reduced  to  less 
than  I  second. 


5.  The  analog  signal  processors  have  large  band  widths  and  very  short  processing  time, 
and  they  serve  as  very  convenient  means  for  detecting  extracted  signal  signatures. 

6.  General  feature  extraction  with  arbitrary  orientation  will  be  possible  if  a  rotational 
capability  of  the  feature  patterns  or  the  sensor  array  is  provided. 


APPENDIX.  PROGRAM  WALSH 


1  FTN4.L 

2  . . .  *WALSH*--REV  1/4/ 78« ••••••*••*•••••• • 

3  C 

4  C 

3  PROGRAM  WALSH 

t  DIMENSION  IA1<  1024),  IA2(  1024  >,  I0ARK!312  ) 

7  DIMENSION  IMPUT!9  ),  LA8IN(  10  ),  LAB0T<  10  ) 

8  DATA  IMPUT/2HDK<2HIN.2HNM<2HPI .  2  HP  T . 2 HT R  >  2HWT /2HTH.2HMC/ 

9  DATA  LAM1N/2NIN.2HPU.2HT  .7*2H  / 

10  DATA  LAB0T/2HTR.2HAN.2HSF  <  2  HO  R  <  2  HR  ,3»2H  / 

11  DATA  1YES/2HYE/ 

12  C 

13  C  GET  LU  OF  FIRST  MCI  AND  OUTPUT  DEVICE. 

14  C  GENERATE  WALSH  FUNCTIONS. 

15  C 

It  CALL  LUNIT!148, IA1. IA2 > 

17  LUIHMAl 

18  HORN*  1 

1*  I T  HR  S  *0 

20  I  X*3  2 

21  IY*32 

22  WR I T  E ! 2  /  3  > 

23  3  FORHATC *LP?*  ) 

24  READ!  I, 30 JICMND 

23  IF  !  ICHHD  E8  .  I  YES  >G0  TO  t 

26  LU0T*2 

27  NCHR  S  *  72 

28  GO  TO  7 

29  6  LUOT-t 

30  NCHRS*1 32 

31  7  1 N  T  V  L  *  1 

32  WRITEC2/8) 

33  8  FORMAT! ’BIPOLAR  INTERVAL?*! 

34  READ!  1. 30 >!CHND 

33  IF  !  1CHND  Efi.lYES )INTVL*-1 

36  CALL  IWTFI  I A 1 .  1 A2, 32, 1NTVL,  32  ,  32, 32  ) 

37  C 

38  C  INPUT  COMMAND  LOOP 

39  C 

40  C  DK--INPUT  ARRAY  DARK  LEVELS 

41  C  IN--INPUT  ARRAY  DATA 

42  C  NM--SET  NORMALIZATION  COHSTANT 

43  C  MC--SET  80UN0S  ON  COEFFICIENTS  COMPUTED 

44  C  TH--INPUT  THRESHOLD  CONSTANT 

43  C  VT--TAKE  THE  WALSH  TRANSFORM 

46  C  PI--PRINT  THE  INPUT  ARRAY 

47  C  PT--PRINT  THE  TRANSFORMED  ARRAY 

48  C  TR--TERH  IN  ATE  PROGRAM 

49  C 

30  C 

31  10  WRITCT2.20) 

32  20  FORMAT !*??•> 

33  READ! 1.30>ICNND 

34  30  FORMAT! A2 ) 

S3  IF  ! 1CNND  EQ . INPUT! 1 >)G0  TO  100 

36  IF  !  !CHND.EQ.INPUT!2))G0  TO  200 

37  IF  !  1CMMD  .  ED  .  INPUT!  3  )  >R£AO!  1 .  •  >NORH 

38  !F  !  ICMND. 64. INPUT!  • DREAD! 1. • >ITNRS 

39  IF  !  1CMHD.E8  INPUT!  9  >>READ!  1.  •  )IX.  IY 


to  IF  <  ICHNO  EO  INPUT! 4  )  1CALL  HA  T  OT  <  I  A  1 ,  L  AB  IN .  3  2  ,  J  2 ,  J  ,  LUO  T  .BCHRS  ) 

tl  IF  ! ICNNO . E8 . INPUT! 3 >>CALL  HO  T  0 1  !  I  A  2  ,  LA  8  Of  .  I  X  ,  I  Y .  3  ,  L  UO  T  .8CH  R  $  ) 

t2  IF  !  ICNNO  E8  INPUT! 7  )>CALL  I V  TF(  IAI<IA2<I2<0<N0RN<  I  X , I Y  J 

63  CALL  PTFOF 

64  IF  <  ICNND  E8 .INPUT!*  MSTOP 

63  GO  TO  10 

66  C 

67  C  INPUT  THE  ARRAY  DARK  LEVELS  AND  PACK  INTO 

68  C  BUFFER  *  1 0 ARK* 

69  C 

70  1  00  CALL  ONAINdAI.  1024. LUIN) 

71  CALL  IPACK( IA1 .  IDARK.  1024  ) 

72  GO  TO  10 

73  C 

74  C  INPUT  ONE  ARRAY  FRAHE,  SUBTRACT  THE  DARK  LEVELS 

73  C 

76  200  CALL  OH  A I  N!  I  A1  ,  1  0  24  ,  LU  1 N  > 

77  J-  1 

73  DO  210  1«1  ,312 

73  CALL  UNPAK! IDARK! I > . IA2, 2 > 

30  IA1C  0  )■  I  A  1  (  J  >- 1 A2 !  1  ) 

31  IF  (  I  A  l  <  J  )  .  L  T  .  0  )  I  A  1  (  J  )  *  0 

82  IF  C  I THRS  EB  0  )G0  TO  201 

83  IF  (IAICJ).LT  ITHRS  )IA1(J>>0 

34  IF  !  I  At  C  J  )  GE  .  I  THRS  )  IA1(  J  )»  100 

33  201  <1*J*1 

36  I  A  1(  J  )*  IA  1  <  J  >-  I  A2C2  ) 

87  IF  (  IA1<J  ).LT.0)IA1!J  )*0 

33  IF  (  ITHRS . EO  0  )G0  TO  210 

33  IF  !  I  Al!  J  >  LT  .  ITHRS  >IA  1!  J  >«0 

30  IF  !  IAICJ  }  .GE .  ITHRS  1IA1C J  >«100 

31  210  J*J*1 

32  GO  TO  10 

33  END 

34  ENDS 


